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Abstract 

We derive an equation for the acceleration of a fluid element in the spherical gravita- 
tional collapse of a bounded compact object made up of an imperfect fluid. We show 
that non-singular as well as singular solutions arise in the collapse of a fluid initially 
at rest and having only a tangential pressure. We obtain an exact solution of Einstein 
equations, in the form of an infinite series, for collapse under tangential pressure with a 
linear equation of state. We show that if a singularity forms in the tangential pressure 
model, the conditions for the singularity to be naked are exactly the same as in the 
model of dust collapse. 



1 Introduction 



The study of gravitational collapse of a bounded spherical object in classical general relativity 
has received major attention over the last few years. The purpose of these investigations has 
been twofold - to establish whether or not naked singularities arise in gravitational collapse 
JTJ, and to study the occurrence of critical phenomena in gravitational collapse 0. 

Various models of spherical collapse have been studied over the last few years, and these 
show that both black holes and naked singularities arise during gravitational collapse. The 
collapsing matter is assumed to satisfy one or more energy conditions, and is in this sense 
regarded as physically reasonable matter. The models studied so far include collapse of dust 
0, null dust 0], perfect || and imperfect fluids ||, and scalar fields 0. In each of these 
cases, the formation of covered as well as naked singularities has been observed. There are 
also demonstrations that both these types of solutions will arise in the collapse of a general 
form of matter ||. Also, a naturalness argument has recently been put forth to suggest that 
covered as well as naked singular solutions arise generically in spherical collapse, subject to 
the assumption of the dominant energy condition [[|. 

Critical behaviour has been discovered in the collapse of matter fields and fluids, largely 
through numerical studies. In some of these studies, it has been observed that the solution 
separating dispersive solutions from collapsing ones is a naked singularity. It should however 
be noted that in studies of critical behaviour, the numerical identification of a "black hole" 
is carried out by the detection of an apparent horizon. The possibility that some of these 
"black hole" solutions are actually naked singularities cannot be a priori ruled out. If the 
singularity is globally naked, there would be a Cauchy horizon lying outside the apparent 
horizon, in the Penrose diagram. Such a singularity cannot be identified in present numerical 
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studies, which do not probe the central high curvature region of the collapsing object. On 
the other hand, a few analytical studies of scalar collapse confirm that the non-dispersive 
solutions contain both black holes and naked singular solutions JTOl . 



In the present paper, we study the spherical gravitational collapse of an imperfect fluid 
under the assumption that the radial pressure is identically zero, but the tangential pressure 
is non-zero. This system has been studied by a few authors in the past [11], and also more 



recently [O, with the recent studies focusing on the issue of naked singularity formation. 



The works mentioned in |T2[ give evidence for naked singularity formation. The purpose of 
the present paper is to demonstrate the occurrence of non-singular and covered and naked 
singular solutions in a specific tangential pressure model. As we will show, our analytical 
results are similar, in some respects, to numerical results of critical behaviour in collapse. 
We consider an equation of state for the tangential pressure pt, of the form px = k(r)p, 
for which we obtain an exact solution of the Einstein equations, in the form of an infinite 
series. For physical reasons, we assume that the function k(r) must vanish at the origin of 
coordinates. As a consequence of this constraint we find that if a singularity forms in this 
model (which is not always the case) the conditions for the occurrence of a naked singularity 
are exactly the same as in the dust collapse model. 

The focus of our investigation will be the singularity that may possibly form during the 
collapse, at the origin r = of the spherical coordinates. This is the so-called central shell- 
focusing singularity. Hence we will specify initial conditions only in a small neighborhood of 
the center and investigate the nature of the collapse in that region, without considering the 
evolution in the other regions of the spherical object. We will also assume that the initial 
conditions are such that shell-crossing singularities do not form during the evolution. 

In Section 2 of the paper we write down the Einstein equations for the collapse of a 
bounded spherical object in an asymptotically flat spacetime. We then derive the equation 
for the acceleration of a fluid element in this spherical object, under the mutual influence of 
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gravity and the two pressures (radial and tangential). In Section 3 we use the acceleration 
equation to derive sufficient conditions for singularity formation, for the cases of collapse 
under tangential pressure, collapse under radial pressure and collapse of a perfect fluid. In 
all these cases we show that there are initial conditions for which the evolution is non- 
singular, and other initial conditions which result in singularity formation. In Section 4 we 
give an exact series solution for the evolution of the area radius, in the tangential pressure 
model under consideration. We show that the same solution also holds for a very specific 
kind of evolution under radial pressure or for a perfect fluid. In Section 5 we give a simplified 
derivation of the earlier results on naked singularity formation in the dust model. In Section 
6 we show that covered as well as naked singularities form in the tangential pressure model 
as well as the radial pressure and perfect fluid models. 



In comoving coordinates (t, r, 9, 0) the spherically symmetric line-element is given by 



where a and u are functions of t and r. The area radius R also depends on both t and r. In 
comoving coordinates the energy-momentum tensor for a spherically symmetric object takes 
the diagonal form T\ = (p, — p r , —pr, —pr)- The quantities p r and px are interpreted to be 
the radial and tangential pressure, respectively. The Einstein field equations for this system 
are 



2 The Acceleration Equation 



ds 2 = e a dt 2 - e^dr 2 - R 2 dtt 2 



(1) 



m! = 4npR 2 R', 



(2) 
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rh = —ATcp r R 2 R, 



(3) 



^=- J ^~ + ^--M-Pr), (4) 
p + p r R{p + p r ) 



. _ 2p AR{p+p T ) 

UJ — — — , \Oj 

P+Pr R{P + Pr) 



and 



m = -R (l +e - a R 2 - e^R' 2 ). (6) 

Here, m(t, r) is a free function arising out of integration of the Einstein equations. Its initial 
value, m(0,r), is interpreted as the mass interior to the coordinate r. 

In order to derive the equation for the acceleration we first define e~ ul R' 2 = 1 + f(t,r). 
Then 



f . 2R' R . 

—uj H = —a 1 (7) 

T3I JDI \ I 



1 + / R> R> 

where the last equality follows from using and eliminating p using (Q) and (^|). Now we 
differentiate (|6|) w.r.t. t after writing it as 



„ • 9 2m 
6 R = li +f 



which gives the acceleration equation 



, m\ , 1, h , e g (l + /) ( 2p' r AR\p T -p r ) \ 

R = -eR(A, Pr + -j + -aR + ^^l-- T -+ . (9) 
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By defining the proper time dr = e a l 2 dt this equation can also be written as 



d 2 R ( m\ (1+f) 



2p' r , 4R'(p T - Pr y 



+ 



(10) 



2R' \ p + Pr ' R(p + Pr, 

For a perfect fluid, we have pt = p r = P, and the above equation becomes (see for instance 
0), 



cPR 
IP 1 



-R Aixp + 



m \ 
R3J 



[1 + f) V' 
R' p + p' 



The Oppenheimer-Volkoff equation for hydrostatic equilibrium is obtained by setting the 
acceleration and the velocity equal to zero, and by noting, from (pi), that / = —2mjR: 



dR 



mp 



i + p - 

P 



1 + 



4nR 3 p 



2m 



(12) 



A few interesting properties about the role of pressure in the acceleration equation (||) should 
be noted. The tangential pressure appears only in the last term, and its gradient does not 
enter the equation. The gradient of only the radial pressure appears in the equation. A 
positive tangential pressure opposes collapse, while a negative tangential pressure supports 
it. 



3 Conditions for singularity formation 

There are various interesting special cases of the Einstein equations for spherical collapse 
given in the previous section, and we consider them one by one. The dust approximation is 
obtained by setting p r = p T = 0. In this case equation (Q) remains as such, while equation 
(0) gives that the mass function does not depend on time t. Equation (f|) implies that a is a 
function only of time; hence we can redefine t and set a = 0. Equation (^) can be integrated 
to get & 1 = R' 2 /(l + /(r)), where f(r) is a function of integration. Hence (0) can be written 
as 

k , = 2n^ +m (13) 

The dust model has been discussed in detail by many authors J|. 
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3.1 Tangential Pressure 



In the case when the radial pressure p r is zero, and the tangential pressure non-zero, consid- 
erable simplification of the full system of equations takes place. As a result of equation (|3]), 
the mass function m(r) is time-independent. Equations @ and (|5D become 

✓ (U, 
The equation (|9]) for acceleration becomes 

+ (16 ) 

We will assume that collapse begins from rest at time t = 0, and choose the initial scaling 
R(0,r) = r, which gives from equation (§) that initially / = —2m(r)/r. In order for collapse 
to begin, the acceleration must be negative, which implies that 

2m 4p T /p 
R l + Ap T /p 

Collapse will continue all the way up to the formation of a singularity R = 0, provided at 
any successive stage in the evolution, the acceleration is negative when R = 0, i.e. provided 
flT7D holds at all later times. The ratio pr/p will in general evolve with time, for a given r. 
Let the initial value of this ratio be denoted by k(r), assumed to be positive. A sufficient 
(though not necessary) condition for continual collapse is that the ratio Pt/p remains the 
same as, or falls below its initial value k(r). If this happens, then ([H]) will be satisfied 
whenever R = 0, because 

2raQ) 2m(r) 4fc(r) > 4p T /p , . 

R r 1 + 4k{r) ~ l+Ap T /p' 1 j 

We now examine collapse with the assumption that the ratio pr/p remains constant during 
evolution, at its initial value k(r). This assumption makes the analysis tractable. The 



constant k(r) is chosen to lie in the range < k(r) < 1. Since the tangential pressure must 
vanish at the origin as a result of isotropy, we must have k(0) = 0. Equation (|]) has the 
solution 



e -(*,r) = 



)p 2 R^ l+k \ (19) 



Here, x{ r ) is an arbitrary function of the coordinate r . Using equation (|2|) we can write 
this as 



e 



-u(t,r) 



xifW R 4k _ R ik _ l + f(t,r) 



Note that C{r) is a positive function. Using the solution (|20D in equation (||) yields the 
following equation for the evolution of the area radius: 



2 

dr J R 

Assuming that collapse begins from rest gives 



dR\ 2m . > t . . 

— = — - 1 + C(r)i? 4fc . (21) 



^(O = (22) 

and hence ( ]21|) becomes 

(f) 2 ^-(-^)(r w 

Let us consider the condition ( |17|) for continual collapse. We are interested in the be- 
haviour near the origin, r = 0. Near the origin, let k(r) behave as a power law, k(r) = ar n . 
Since 2m /r goes as r 2 near the origin, collapse cannot take place if n = 1, and will neces- 
sarily take place if n > 3. The case n = 2 is critical. Now, collapse will take place provided 
a < 27rpo/3, but not otherwise. Here, po is the initial central density. Thus we find that for 
positive k there are singular as well as non-singular solutions. In Section 5 we will discuss 
conditions for the singularities to be covered or naked. 
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3.2 Perfect Fluid 



The perfect fluid approximation is obtained by setting p r = p T = p. Equations @, (Q) and 
(|6|) remain as such, while equations (f|) and become 

✓ = (24) 

P+P 

2 P 

u = -. 25 

P + P R 

The acceleration equation becomes 

* R -- B (^ + ™\-l±+llj!-. (26) 



dr 2 V i? 3 / i2' p + p' 

For collapse starting from rest at time t — and by using the initial scaling R = r we get 
that initially 1 + / = 1 — 2m(r)/r. In order for the acceleration to be negative initially we 
need 

2m > 1 + 4?rpr(p + p)/p / 
r 1 — (p + p)/2rp' 

Let the initial density profile near the center be 

p = p - ar 11 , a > 0, n > 2 (28) 

and let the equation of state be p = kp, with k a positive constant. For n = 2 the condition 
for negative initial acceleration is 

a < 27r(k+ l)p (l + Po/3fc)- ( 29 ) 



For n > 3 the initial acceleration is necessarily negative. It can be checked that a sufficient 
condition for singularity formation is that the quantity pR 2 increases over, or remains the 
same as its initial value, and further, the quantity R\dp/dR\/(p + p) decreases from, or 
remains the same as, its initial value. 
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3.3 Radial Pressure 



Corresponding equations may be written for the case when the tangential pressure pt is 
zero, and collapse takes place only under radial pressure. It can be shown that starting from 
rest, the fluid will undergo collapse leading to singularity formation, provided the following 
condition is satisfied initially and subsequently: 

2m 87T Pr R 2 + [2p' r R + 4p r R'} / (p + Pr )R' 

R > -l + [2p' r R + 4p r R>] /{p + Pr )R' ' { } 

Consider the initial density profile near the center to be of the form (|28|) (now with n > 1), 
and an equation of state p r = k(r)p. Since the radial pressure must vanish at r = 0, k(r) 
must vanish at the origin. Let k(r) = Ar n near the center. Then it follows from equation 



(|30|) that collapse will not take place for n — 1, will necessarily take place for n > 3, and 
if n = 2, collapse will take place provided po > tt/QA. A sufficient condition for singularity 
formation is that the quantity on the right hand side of ( |3~0D remains constant or falls below 
its initial value. 



4 An exact solution 



We will be interested in solving the collapse equation (|23"D, which describes the collapse of a 
fluid starting from rest, subject to a tangential pressure equation of state px = k(r)p, with 
k(0) = 0. If this equation can be solved for a given m(r), the remaining unknown functions 
(p, u), a and px) can also be obtained. 

Interestingly enough, there is also a class of evolutions for the perfect fluid case and the 
radial pressure case, for which the equation for evolution of the area radius can be cast in a 
form similar to (^). Consider evolutions of the kind 

R' 

8irpR 2 = 6(r), a' = ^(r)—. (31) 

R 
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One cannot attach any physical importance to these assumptions, and they may at best 
only be approximately obeyed during the evolution. Their advantage is that one can then 
obtain a solution, subject to these assumptions, and demonstrate critical behaviour, and the 
occurrence of covered as well as naked singularities, in this solution. Our main purpose is 
to give a solution for the tangential pressure problem. It so happens that on the side, we 
can say something useful about a specific perfect fluid model, whose evolution equations are 
similar to the tangential pressure model. However, the fact that an exact solution is possible 
could be a motivation for understanding the assumptions on physical grounds. We would 
like to point out the following in that context. 

It can be easily seen from equation (|I|) that a plays the role of the Newtonian potential 
in the weak field limit, suggesting that a' represents the gravitational pull on the source 
matter. One could imagine a fluid source for the Einstein field equations wherein the fluid 
particles were moving as shells of constant r. The mean curvature of these shells would then 
be R' /R. Also, 87tpR 2 is analogous to the force experienced by the body of fluid enclosed by 
the shell, solely due to the dynamics of the source particles. 

With the assumptions made in (RT]) we can integrate equation (S) to get 



where m (r) is the initial mass distribution. Similarly, (0) and (|3lD can be used to get 



where A(r) is an integration constant. Assuming that collapse begins from rest at t — 
(where the scaling is R = r), we get that A(r) = (1 — 2m (r)/r) r~^. 




(32) 




(33) 



To rewrite equation (|6|), the following transformation is made use of: 




(34) 
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where the requirement of exactness of the differential equation restricts the choice of Z(r) 
to 



,<t' 



Z{r)= Je^- + g(r) (35) 

g being an arbitrary function of r. We aim at solving equation (§) keeping r fixed. We vary 
t under this restriction and examining equation (|B|) obtain 



2 

dR\ 2mo / ' r \ ( 2mo(r 



-i + e(r) --i -P^-i 3 ■ (36) 



(If 



In case of a perfect fluid, the assumptions in ([31]) reduce to 

8,pR^ = e(r),-^^ = ^r). (37) 

For an equation of state p = kp and an initial density profile given by near the 
center, we get that near r = 0, 

, / % nak 

*W ~ (ITT)ft (38) 

After solving (|36"D for i?(r, r) we can obtain the solutions for (p, u, a and p#). 

For collapse under radial pressure the assumptions fl3~T|) reduce to 

2 + gp^ 

STrpi? = 0(r), R'(p + p ) = ^( r )- ( 39 ) 

For an equation of state of the form p r = k(r)p, with k(r) = Aor n , we get that ip(r) = 
-4(n + 2)k(r). 



The equation (|23|) for evolution of the area radius in the case of collapse with tangential 
pressure is a special case of (|36|) , obtained by setting 6{r) = and if)(r) = 4k(r). We now 
obtain an exact series solution of Equation (j36|), for a general 9(r) and ip(r). After defining 
the variable y = R/r we write ( |36|) as an integral: 

dr = - / f 1 - - 1+6 , y + 1 ~ 2m °/y +^ ~ 1/2 . (40) 

' '0 + 2m /r V + 2m o /r y # + 2m /r ^ V ; 
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By using the binomial expansion 



. 1/2 ~ 1.3.5...(2n-l) 
(1 - X) 1 = 1 + -ar 



n=l Z 



and by taking 



1 + g + 2m /r - 1 
+ 2m /r 6 1 + 2m / r 



the integral becomes 

f*r = - , r \ 2 -y^ + £ L3 - 5 - (2 r 1) f^Vvdy 



s /9 + 2m /r 
We write x n as 



and hence get the solution to the integral as 



^ (2n)\ 1 



2m \ j y n +^+ 3 / 2 



£ n Cj (1 + 6)"- j — - - 1 



r / n + jvp + 3/2 



3=0 

This is an exact solution of the Einstein equation (pq) . 

We cast this into the following form for convenience, 

r - ro(r) = - r y 3/2 G(2/,r) 
^0 + 2m /r 

where 



^ o 22n( n )2(fl + 2mo /r)™ X 



2m V y n +M+ 3 / 2 



V n C ? (1 + 0) n ~ J - 1 , 



13 



It is useful to work out the derivatives of G : 

3G(y,r) I 



dyJr 



y 



+ 



2m /r + 9-(l + 9)y- (2m - l)y^ +l 
2mo/r + 9 



(49) 



'dcr 

dr 



E 



(2n)\ 



X 



a ^ 2 2n («0 2 (0 + 2m /r)' 



j=0 



(50) 



where 



. 0' 2mUr — 2m /r 2 2mUr — 2m /r 2 + 6' 
D{y,r,n,j) = (n - j)^—^ + j — ; n + 



1 + 



2m /r - 1 2m /r + 6 

jif)' \og(y) - 



(51) 

n + jip + 3/2' v ; 

For the special case of dust, we have 6 = ip = 0, and r — t. It can easily be shown from the 
dust equation ( |T3"D for the case of collapse starting from rest that the solution is given by 

R 3 / 2 farcsin 



t - t (r) 



^2m(r) 



y 



3/2 



V 



(52) 



The series solution fl4"S ) can be summed up in the dust case and it can be shown that the 
sum is equal to the closed form solution given in fl52|). 



5 The Dust Solution 



We will be interested in using the solution (|46| ) to find out whether the singularity that 
forms at r = is covered or naked. We start by noting that in the cases of interest discussed 
above, ip(r) has a power law form, near r = 0. Also, 0(r) has a power law form near r = 0. 
Because of this, it can be inferred from ( fl6|) that as r — > the solution approaches the dust 
solution given in (|52|) , with the difference that in (0), t should be replaced by r and 2m(r) 
should be replaced by 2mo(r) + r0(r), during the approach. 
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Hence it should be possible, as explicitly shown later, to use the naked singularity analysis 
for dust, carried out earlier, to draw conclusions about the occurrence of a naked singularity 
in collapse of a fluid under tangential pressure or radial pressure, or a perfect fluid, subject to 
the equation of state chosen above. Before we do so, we would like to use this opportunity to 
present a simplified derivation of the dust naked singularity. We give this simplified analysis 
first for the marginally bound case f(r) = 0, and then for the general case. 

The solution to equation ( |T5D in the marginally bound case is 

^3/2 =r 3/2_3y^-^ (53) 

An initial scaling R — r at the starting epoch t = of the collapse has been assumed. A 
curvature singularity forms at r = at time to = 2/3^/po, where po is the initial central 
density. We assume a series expansion near r = for the initial density function p(r), 

p( r ) = Po + Pir + 7yp2r 2 + ^p 3 r 3 + ... (54) 

The series expansion for the mass m(r) can then be deduced using (||). 

From equation ([53), we evaluate R! and then substitute for t from the same equation. 
In the resulting expression, substitute R = Xr a , and perform a Taylor expansion of F(r) 
around r = 0, so as to retain the leading non-vanishing term. We then get that near r = 0, 



JL=X+ ^ r «H-3/2-3a/3_ (55) 



Here, (3 = —2qm q /3mo, and q is defined such that in a series expansion of the initial density 
p(r) near the center, the first non- vanishing derivative is the qth one (=p q ), and m q = 
4:7ip q /(q + 3)q\. It may be shown that a has a unique value at the approach to the central 
singularity, given by setting the power of r to zero in the second term, i.e. a = 1 + 2q/3. 
This reproduces the result of the R' calculation performed earlier [H|, in a simpler manner. 
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For the non-marginally bound case (i.e. / ^ 0) the solution of the Tolman-Bondi equation 

is 



R 3/2 G(-fR/F) = r 3/2 G{-fr/F) - J¥{r)t 



(56) 



where F(r) = 2m(r), and G(y) is a positive function having the domain 1 > y > — oo and is 
given by. R' can be evaluated as before, and then we eliminate t and substitute R = Xr a . 



The power-series expansion for f(r) near r = is of the form [14 



f(r) = f 2 r 2 + f 3 r 3 + Ur 4 + ... (57) 
which implies that the argument (—fR/F) of G on the left-side of fl56|) goes to zero as r — > 0. 



The derivative of G(—p) w.r.t. its argument p = rf/F is obtained by differentiating ( |56|) 
which gives 



dG(-p) 3G 



d(—p) 2p py/1 +p 
Using this, one can now perform a Taylor-expansion of G, F and / about r = to get exactly 
the same expression for R' as given in (|55|) above, where (3 is now given by 



P = q 1- 



A 
2F n 



G{-h/F ) 



F 



2/2 



1 + 



A 

2F n 



+ 



q+2 



(5f 



The constant g is now defined such that the first non-vanishing derivative of the initial 
density is the gth one, and/or the first non- vanishing term in the expansion for f(r), beyond 
the quadratic term, is of order r q+2 . The constant a is again equal to (1 + 2q/3). Thus the 
R' calculation is simplified for the non-marginal well. 

In order for the singularity at r = to be naked, radial null geodesies should be able 
to propagate outwards, starting from the singularity. A necessary and sufficient condition 
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for this to happen is that the area radius R increase along an outgoing geodesic, because R 
becomes negative in the unphysical region. Thus we write, along the geodesic, 

R N 



dR 

dr 



R' + R^ 

dr 



R' 1 + 



R! 1 



' f + F/R 
l + f 



(59) 



Here we have substituted dt/dr = R' /\/l + / along an outgoing null ray (using the metric 
(|l|)) and substituted for R from (|13|). dR/dr should be positive along the outgoing geodesic. 
We now define u = r a , and use X as a variable, instead of R. Hence, in the approach to the 
singularity, (i.e. as R — > 0, r — > 0), X takes the limiting value X given by 



lim — 

r->0,fl-«-0 u 



lim 



dR 



ar 



a-l 



dr 



lim 



ar 



a-l 



R! 1 



I f + F/R 
l + f 



(60) 



By using (55) we can write 



a \ x/X 




/(0)+A /X 



(61) 



\ l + /(0) 

The constant A is the limiting value of A(r) = F(r)/r a as r — ► 0. The variable X can 
be interpreted as the tangent to the outgoing geodesic, in the R, u plane. As can be seen 
from the above, the positivity of dR/dr along an outgoing geodesic is equivalent to requiring 
that the equation ([H]) admit a positive root X . This will depend on the initial density and 
velocity distribution, which determine the functions F(r) and f(r), and hence the functions 
(3 and A. One can solve equation (pip, and as shown in |14|, the results for the nature of the 
singularity are the following: 

In the marginally bound case, the singularity is naked if p\ < 0, or if p\ = 0,p 2 < 0. 
If Pi = and P2 = 0, one defines the quantity £ = 2m 3 /(2mo) 5 ^ 2 . The singularity is naked 
if £ < —25.9904 and covered if ( exceeds this value. If p\ = p 2 = p 3 = 0, the singularity 
is covered. In the non-marginally bound case, if F\ and f% are non-zero, the singularity is 
naked if 



Qi 



l 



A 

2F n 



G(-f 2 /F ) 



Fi_3fc 
F 2/ 2 



1 + 



A 

2F n 



+ 



FY 
Fo, 



(62) 
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is positive. If F\ and / 3 are both zero, and F 2 and f^, are non-zero the singularity is naked 
if 



Q2 



1 - 



2F , 



G(-f 2 /F ) 



1 + 



/a 



2F 



fi F Q 



(63) 



F2 _ 3/4 
F 2f 2 _ 

is positive. If F%, / 3 , F 2 and / 4 are zero and if F 3 and / 5 are non-zero, the singularity is 
naked if 



1 



A 

2F n 



G(-f 2 /F ) 



F3_3fc 
F 2f 2 



1 + 



A 

2Fn 



_|_ /jj _ -^3 
/2 ^Q, 



(64) 



is positive. If Fi, / 3 , F 2 , / 4 , F 3 and ^5 are all zero, the singularity is covered. 



A special case of non- marginally bound collapse is the collapse starting from rest, for 
which f(r) = —2m(r)/r. As a result, G(—fR/F) = tt/2 and the calculation of R' in (|56|) 
can be carried out exactly as in the marginally bound case, to get the result (|55|) , with 
(3 = — Tiqm q / Auiq. The conditions for a naked singularity to occur are the same as those 
stated above for marginally bound collapse, except that £ is defined as £ = 67rm 3 /4(2mo) 5//2 . 



We also point out that in a recent work JOJ we have given a yet simpler derivation of 
the dust naked singularity, which directly looks for a self-consistent solution of the geodesic 
equation, in a neighborhood of the singularity. In principle, the method described in jOJ 
can be applied also to the tangential pressure solution discussed here. 



6 The occurrence of covered and naked singularities 



We can now utilize the dust results for inferring the occurrence of covered and naked sin- 
gularities in collapse with the equations of state considered here, except that two further 



subtleties remain to be sorted out. Firstly, the solution ([46]) is written with respect to the 
time variable r, for a fixed r. We need to show that a relation similar to equation ( |6T)D holds. 
Secondly, for a general ip and 9 (i.e. not restricted to power law forms for instance) the 
dependence of / and F on t as well introduces terms in X different from the dust case. It 
is not obvious a priori that these vanish as one approaches the singularity, even though the 



solution tends to behave like dust in this limit. It needs to be explicitly shown that a rela- 
tionship identical in form to equation fl6l] ) holds in the general case. This is demonstrated 
in the rest of the section. 



By requiring that dr be integrable, we can write 



dr = e a/2 dt + Z(r)dr = e a/2 dt + 



dr 



(65) 



where the last expression follows from using the integrability condition on Z(r), and g(r) 
is an arbitrary function. We choose g(r) in such a way that in the approach to the central 
singularity the coefficient of dr in (|65|) vanishes. As a result, in the approach to the singularity 
we can write dr = e a l 2 dt. 



The following can be easily shown from ( p6|) and the solution 



'dR\ 



l -[A{y,r)] 1/2 



where 



dR/dr is finite in the approach to the singularity. Also, 



(66) 



(67) 



dFC 
dr 



1 + 



2mo jr + 9 



yjA(y,r)B(y,r) 



where 



B(y,r) 



t qG _G + 2m' /r-2m /r 2 + 6' _ fdG\ 
t — tq r ' 2iriQ jr + 6 \dr J 



Now, 



(R') t = (R% + 



'dR\ (df 
dr \ dr 



(68) 



(69) 



(70) 



and it follows from the above discussion that the last term, M(y, r) say, can be made to 
vanish in the approach to the singularity, so that in the limit the quantities (R')t and {R') T 
become identical to each other. 
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From assumptions ( |31"D one can easily show that 



R l e -" = (1 - 2m /r) y* 



(71) 



This leads to 



j/2 y_ 



l-V/2 



1 — 2mo/r 



1 + 



2mo /r + 9 



\l A{y,r)B{y,r) 



.l-V/2 



+ 



1 — 2mo/r 



M{y,r) (72) 



From fl6"5|) we also get that along a null geodesic 

% = e^ + Je^ a -dt + g (r). 
This, after using (72), is re- written as 



9r - 2m /r 



^A(y,r)B(y,r) 



+ 



V 



1-^/2 



(73) 



Y^l - 2m /r 



M(y,r) 



(74) 



Hence we may examine the rate of change of i? along an outcoming null ray as 



v / riuZZ geodesic \ / 



nuZZ geodesic 



'dr 
dr 



(75) 



From equations (|6q) , (|68|) and (|74l), this can be cast as 
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(76) 
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where the derivative on the left is to be evaluated along null geodesies. This resembles the 
roots equation (|59|). 

Indeed, the first term with square parenthesis in equation (|76|) is similar to R! jr a ~ x 
in equation fl60l) . It turns out that as one takes the limit r — > 0, the third contribution 
within this term containing the summation vanishes for physically reasonable choices of ijj, 
i.e. vanishes as one takes the limit. This requirement is certainly satisfied by the three 
special cases in Sections 3.1, 3.2 and 3.3. In the rest of the contribution, these conditions 
imply that the result is the same as dust, save the replacement of 2m by 2mo + r9 and to 
by To in the approach to the singularity. Similarly, the second term in square parenthesis 
in equation (|76| ) has a counterpart in the dust case, although the familiar / and F of (|59| ) 
have no individual correspondences. However, if one examines the limiting form fl60|) of dust 
and works out flTS| ) in the limit, then one finds that /(0) of ( pC| ) corresponds to — 2mo/r — 9 
evaluated in the limit r = and A gets replaced by (2m + r6)/r a (evaluated in the limit 
r = 0). Hence we arrive at the same roots equation as ( ftTl) . 

It is important to realize that the first subtlety mentioned at the beginning is sorted out 
by a judicious choice of the r coordinate and results in the fact that (dR/dr) T tends to R! of 
the dust case in (f>ty) in the limit. The second subtlety presents itself explicitly when it turns 
out that the correspondences of / and F of dust are not f(r,t) and F(r,t) of the general 
case if one has not taken the limit. It would be therefore in general incorrect to conclude 
equation (f76|) to be equation (|59f) with / and F simply replaced by their generalizations, at 



this stage. However, the term (/ + F/R)/(l + f) in fl5§D still corresponds to the quantity 
obtained by simple generalization of the free functions of dust. The existence of such a 
quantity is not apriori guaranteed when one notices that the solution (|46|) behaves similar 
to the dust solution in the approach to the singularity. 

We now discuss the nature of the singularities. Consider first collapse under tangential 
pressure pr with the equation of state pr = k(r)p, which we have discussed above. Since in 
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this case 9{r) =0 it follows from the exact solution given above in ( f4"6| ) that if the collapse 
ends in a singularity, then the conditions for the occurrence of covered and naked singularities 
are exactly the same as in the case of dust collapse starting from rest. The introduction of 
a tangential pressure (which must vanish at the center even though it is non-zero elsewhere) 
does not change the nature of the dust singularity. Further, we expect that even if we do not 
restrict to the case of linear equation of state, the conditions for a singularity to be naked 
in the tangential pressure model will be exactly the same as in the dust case, so long as the 
tangential pressure vanishes at the origin. 

The situation is more interesting in the perfect fluid case, because now 9(r) ^ 0. The 
leading order solution is the dust solution, but with an "effective mass" 2mo(r) +9(r). Hence 
the dust results on the nature of the singularity are applicable, except that we must replace 
2mo(r) by 2mo(V) + 9{r). Recalling the definition of 9{r), we can expand it in a series near 
r = 0, 

°( r ) 2 1 2 1 6 

47T 2 o 

We have also assumed an equation of state p = kp, with k a constant. Hence the coefficients 
in the above expansion for the pressure are related to those for the density in the expansion 
in ([54]). Eqns. (^) and (|29| ) give the conditions for singularity formation, if p 2 < 0. If a 
singularity does form it will be naked. If p2 = and ps < a singularity will necessarily 
form. It will be naked if ( = 37r(2m3 + 27rp3/3)/4F 5//2 is less than —25.9904 and covered if 
( exceeds this value. Thus we find that in the case of a perfect fluid, the condition for the 
occurrence of a naked singularity differs from that in the dust case, because of the presence 
of the constant k in the definition for By following a similar series of arguments one 
can conclude that if a singularity forms in the radial pressure model considered above, the 
conditions for it to be naked or covered are exactly the same as in the dust case. 
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7 Appendix 



Although we have shown that the collapse of the cloud in Sections 3.1, 3.2 and 3.3 leads 
to the formation of a singularity for some initial conditions, it is necessary to ensure that 
the same happens in case of the general solution ( [46]) obtained with the assumptions (3T[). 
We proceed in a manner similar to these cases. Initially, let R be and R be negative. If 
the collapse has to take place without any rebounds hereafter, one needs to ensure that the 
shells have negative 'acceleration' R whenever they reach the velocity limit R — 0. This 
along with the initial condition keeps R non-positive. It is straightforward to show from 
assumptions ([H]) using ( |71~D and @ that this requirement for avoiding rebounds reduces to 
the following condition on R: 

R < f (2m /r + 9) 



r (2m /r + 9) + (1 + 9) f 



whenever R satisfies 



2m /r + 9- (1 + 0)— + (1 -2m /r) (-) =0. 

Note that this is equivalent to an algebraic inequality dependent solely on the initial data, 
when the root R/r in the equation above is bounded by the first inequality. Let us suppose 
this constraint on the initial free functions is satisfied for one of the roots which we know to be 
1. This is nothing but the initial epoch of the collapse. The collapse begins. Let us suppose 
there occurs another real root (i.e. R vanishes again). This time the collapse has proceeded 
and the root, therefore, has to be lesser than 1. The first inequality, already satisfied for 
R/r = 1 will therefore be automatically satisfied making the acceleration negative in this 
situation. Thus the rebound will be prevented. Hence we conclude that if we ensure that 
the initial data is chosen such that the cloud begins to collapse at the first instant, then it 
is implied that the cloud will have no rebounds at all at any later instant. 
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